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We study deformation of a cavity around a large projectile moving with subsonic 
velocity in the cloud of small dust particles. To solve this problem, we employ the 
Navier-Stokes equation for a compressible fluid with due regard for friction between 
dust particles and atoms of neutral gas. The solutions shows that due to friction, the 
pressure of dust cloud at the boundary of the cavity behind the projectile can become 
negative, which entails formation of a microscopic void free from dust particles - the 
cavity deformation. The corresponding threshold velocity is calculated, which is 
found to decrease with increasing cavity size. Measurement of such velocity makes 
it possible to estimate the static pressure inside the dist cloud. 

PACS numbers: 52.27.Lw, 83.10.Rs, 82.70.Dd 
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I. INTRODUCTION 



Complex (dusty) plasma is a low-temperature plasma including microparticles. Due to 
electron absorption, particles acquire a considerable electric charge. Thus, a strongly coupled 



Coulomb system is formed 



49]. Such plasma represents a natural system which makes 
it possible to observe various collective phenomena at the level of individual particles. In 
experimental setups, complex plasmas are usually studied in gas discharges at low pressures, 
e.g., in radio frequency (RF) discharges. A large homogeneous bulk of complex plasma, 
which almost fills the entire discharge volume, can be observed under microgravity conditions 
either in parabolic flights 
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19|]. 



14j or onboard the International Space Station (ISS) 
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In recent studies, attention was focused on tracer particles or projectiles moving through 



a cloud of complex p 



of acceleration 
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asma. Such projectiles are generated using controlled mechanisms 
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2l|. Projectiles moving with 



20j . or they can appear sporadically 
supersonic velocity lead to the formation of extended Mach cones, subsonic (slow) ones 
produce localized disturbances of surrounding particles. In Ref. Q, it was suggested that 
the latter regime, realized when a relatively large subsonic projectile moves in a dense cloud 
of smaller particles, can be well approximated by a flow of an incompressible fluid. In the 
framework of the same hydrodynamic approximation, it was demonstrated that, along with 
the regular neutral gas drag, there is an additional force exerted on the projectile due to 



friction between neutral atoms and the particle fluid 



23) 



In this study, we try to estimate the deformation threshold for a (initially spherical) cavity 
around a subsonic projectile. This deformation e merg es abruptly as the projectile velocity 



increases, as it is seen in snapshots given in 
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To solve this problem, we employ 
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the Navier-Stokes equation for a compressible fluid taking into account friction between 
dust particles and atoms of neutral gas. We do not treat the deformation self-consistently. 
Instead, we imply that the projectile velocity is below the threshold value, so that a regular 
flow around a spherical cavity with no stall can be treated. The solutions shows that due to 
friction, the pressure of dust cloud at the boundary of the cavity behind the projectile can 
become negative, which entails the formation of a microscopic void free from dust particles 
- the cavity deformation. This occurs at some threshold velocity which decreases with 
increasing cavity size. Measurement of such velocity would make it possible to estimate the 
static pressure inside the dust cloud. 

The paper is organized as follows. In Sec. [TTJ, we solve the non- stationary Navier-Stokes 
equation for an incompressible particle fluid, in Sec. IHIl the gas dynamics problem is solved 
and the corrections for fluid compressibility to the velocity and pressure fields are calculated. 
Calculation results are compared with available experimental data in Sec. IIVI the results of 
this study are summarized in Sec. [V] 



II. INCOMPRESSIBLE FLUID APPROXIMATION 

Consider an irrotational flow of incompressible particle fluid formed by the dust crystal 
melted around a projectile moving with the velocity u(t) relative to the dust, where t is the 
time. In this Section, we will assume that the cavity around a projectile is not deformed 
and has the form of a sphere with radius R. Note that this is the only possible form of a 



cavity when the parameter of interaction between the projectile and dust particles 



is great (under typical experimental conditions, (3d p ~ 100 



22|, 



22j). The Navier-Stokes equation 



describing the velocity field v(r, t) in the reference frame of the projectile has the form 



— + V- V)v + z/(v + u) = ^ 

ot p 



(1) 



where p = M^n^ is the mass density of dust fluid (assumed in this section to be constant), 
with M d and being, respectively, the mass and number density of dust particles of the 
radius a<i- Furthermore, p(t, r) is the dust pressure field, v = (8\/2Tr /3)5m r ,n v .VTn.a^ / Mh. 



24j, and m n , 



is the friction coefficient with 5 ~ 1.4 being the accommodation coefficient 
n n , T n , and vt h = (Tn/mn) 1 ^ 2 are the mass, number density, temperature, and thermal 

rom the 



velocity of neutral gas molecules, respectively. Equation pj) assumes also that far 
projectile, dust particles are quiescent relative to neutral gas. It was shown in |22j that 
the approximation of nonviscous flow results in a fairly good description of trajectories of 
individual dust particles; estimate of the viscosity term for the dust particle fluid is indicative 



of the fact that in most cases, it is small 23]. This allowed us to omit it in Eq. ([T]) and to 
confine ourselves to the nonviscous approximation. 

For an incompressible fluid, the continuity equation is reduced to 



V- v = 0. 

The boundary conditions for ([1]) and (j2J) are (v ■ n)| r=JJ = 0, where n = r/r, and v 
at r = oo. For an irrotational flow (V x v = 0), we substitute v = Vy? — u in Eq. 
obtain the equation 



VV = 0, <p{t, oo) = 0, 



which has the solution 



25) 



dip 
dn 



0. 



--R 



R 3 



or 



R 



v r 



2r 3 



- [u — 3n(n • u)] — u. 



(2) 



u 



to 



(3) 



(4) 



(5) 
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The pressure field p(r) can then be found by substituting velocity in Eq. ((T]). This 
yields the following pressure distribution at the spherical surface of a cavity (r = R), 



where p st = const > is the static pressure of dust at r = oo. 

Consider deformation of a cavity around a projectile propagating in the dust crystal 
(Fig. [2J. Here, we will neglect finite fluid compressibility - corresponding corrections are 
discussed in Sec. IIHI Deformation can be caused solely by formation of a void in the stream of 
dust fluid because the interaction parameter /3^ p is large. Since no evidence for an appreciable 
attraction between dust particles has been reported in the literature, we model the cloud 
of dust particles by a system of repulsing soft spheres. In such a system, a cavity can be 
deformed by void formation if the condition p < is satisfied over some area around a cavity. 
This condition corresponds to the cavitation condition in a metastable fluid and to stall in 
gas dynamics. One can associate the shape of a cavity with the line p(r) = 0. 

Let us first consider the case u — 0. According to Eq. (jSJ), friction between dust particles 
and neutral gas gives rise to pressure increase in front of the projectile and its decrease 
behind it, so that p(r) reaches a maximum at the "front pole" (where u • n = u). By 
introducing the cosine of the polar angle, cos# = u ■ n/u = £, we readily derive 



where £ cr = — — (uR/u). If |£ C r| > 1 then the pressure reaches a minimum at the "rear 

1 9 

pole" (£ = —1), where p min = p st — -pw 2 (-|£ cr | — 1)- This corresponds to the low-velocity 

2 

regime u < g u Ri when the minimum pressure decreases (approximately) linearly both 
with u and R. The cavity deformation threshold is defined by the condition p m i n = 
and is reached easier for larger projectiles, at the critical velocity u cr ~ 2p st /{ppR) [here 






(7) 
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1 2 
we assume p st < -p(uR) 2 ]. Note that in the high- velocity regime u > —vR the pressure 
8 9 

minimum is shifted to a certain latitude in the rear hemisphere, approaching the "equator" 
(£ = 0) asymptotically. In this case the friction plays minor role and the critical velocity is 

u CT 



It is well known in hydrodynamics that the presence of points or small regions of negative 
pressure can lead to cavitation. However, in complex plasmas the distance between fluid 
particles is likely to be even greater than the size of such regions. Therefore, the continuity 
approximation fails to describe the cavitation in complex plasmas and, hence, this process 
is unlikely to be observed in a real experiment. On the other hand, upon further increase 
of the projectile velocity the regions of negative pressure should grow until they cover an 
appreciable part of the surface of a cavity, including its rear pole. This should be followed by 
a stall observed as cavity deformation. At this stage, a growing deviation from the spherical 
geometry makes analytical solution of the problem too complicated. Moreover, to our best 
knowledge, there is no criterion that would tell us how large the negative region should be to 
initiate a stall. Nevertheless, one can speculate that the velocity u cr derived above provides 
a reasonable order-of-magnitude estimate for the cavity deformation onset, so that one can 
also estimate the static pressure p st . 



Note that for arbitrary acceleration (assuming u||u), one can introduce the effective 
damping rate z/ eff = v + ii/u which can have either sign, so that the critical cosine £ cr can 
have either sign as well. Thus, if a projectile experiences sufficiently strong deceleration 
{\ii\/u > v) the cavity can be also deformed at the front hemisphere. 
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III. FLOW WITH FINITE COMPRESSIBILITY 

In this Section, we will estimate the effect of finite compressibility of the dust cloud 
on the velocity and pressure fields. Here, we will treat a steady (u = 0) and irrotational 
(V x v = 0) flow, so that the fluid density is a function solely of the radius- vector: p = p(r). 



Since V x v = and Vp = c 2 Vp, where c = J (dp/dp)T is the sound velocity, we can rewrite 
Eqs. ([I]) and (jSJ) in the form 



c 2 



(v • V)v + i/(v + u) = -— Vp. (8) 



In this case, the continuity equation includes the fluid density, 



V-v + v^ = 0. (9) 



Here, c is assumed to be constant independent of the local velocity and plasma state param- 
eters. We substitute Vp/p = — (l/c 2 )(v • V)v — (z//c 2 )(v + u) from (jSJ) in (jSJ) to derive 

Vv = -^v (v V)v + -^(v 2 + u- v). (10) 

We will search for the solution of Eq. ffTU]) in the form v = v + v 1; where v (r) is the 
solution of incompressible problem (for p = p = const, i.e., for c — >■ oo) given by Eq. (j5J), 
the respective pressure distribution po(r) is presented by Eq. (jfJJ) (for u = 0). Obviously, 
V • vq = 0. We assume that |vi| |vq| and retain solely zeroth-order terms on the rhs of 



(flQj) to obtain 



V-vi = ^v - Vv 2 + ^(v 2 + u-v ). (11) 



Here, we took into account that V x v = and hence (v • V)v = Vf 2 /2. Since |v | ~ u, 
|V| ~ 1/R, and {u/c) 2 is small for subsonic motion, the estimate |vi| ~ (u 2 /c 2 ) max{«, uR} 
following from ffTTl) justifies our assumption. It is seen from this estimate that both sides 



8 

of Eq. (flTj) are of the order (u/c) 2 , and terms other that zeroth-order ones would lead to 
higher-order terms in this parameter. Obviously, |v x | — > at c — > oo. 

We introduce the scaled coordinate x = r/R and the dimensionless potential (p\ defined 
by the relation Vi = (u 3 /c 2 )V<^i (in so doing, the condition V x V] = is satisfied) to 
rewrite (fTTl) in the dimensionless form 

-Ki"^^-^^^)^ < 12) 

where -P;(£) are the Legendre polynomials and k = vR/u. Here and in what follows, dif- 
ferentiation with respect to scaled coordinates is assumed. The boundary conditions for 
Eq. (fl2l) follow from the conditions Vi(oo) = and n ■ vi | r =R = 0, i.e., 

dip! 



<£i(oo) = 0, 



dn 



0. (13) 



x=l 



Solution of Eq. f[T2"|) is derived in Appendix |A] The resulting finite-compressibility cor- 
rections to the potential and velocity at the cavity surface follow from ( 1A6I) : 

81 3 k /l 2 1\ 551 
<Pi \x=i - y^£ + 2 U ~ 3/ ~ 880' ( U ) 

u 2 /243 2 551\ , , V ; 

Vi L_i = — £ H — £ (u — fun). 

1 |!B-1 c 2 V176^ 4^ 880/ V ^ ; 

As it could be expected from symmetry considerations, vi \ x =i = at £ = ±1, i.e, at both 
poles. Note that result ( 1T4"|) validates the estimate |vi| ~ (u 2 /c 2 ) max{«, zaR} given above. 

On the basis of velocity correction ffT4"]) and using Eq. (jSJ), it is straightforward to calculate 
a compressibility correction to the pressure, p\ = p — po- This is performed in Appendix 
IB1 It is worth mentioning that in the general case, correction to the pressure (IB4[) includes 
solely even powers of £, so that corresponding correction to an additional force exerted 



on the projectile due to dissipation in the surrounding particle fluid [Eq. (10) in Ref. 



23|] 
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vanishes. In the friction-dominated regime (large k), one can retain in (IB4j) solely the term 
proportional to u 2 , which is independent of £, so that the total pressure is 

PoUvR p u 2 ( a 4zAR2\ 
P(0 - Pst + — g— i + — I 9f - 5 + I . (15) 

Thus, the pressure increases due to finite compressibility, so that the threshold velocity of 
cavity deformation increases as well. On the other hand, one can see that the resulting 
correction is small when k(u/c) 2 < 1. As we show in the next section, the compressibil- 
ity correction can usually be neglected for subsonic projectiles under typical experimental 
conditions. 



IV. COMPARISON WITH EXPERIMENT 



Based on the developed theory, let us analyze the possibility of cavity deformation in 
typical experiments with complex plasmas. We first focus on the experiments carried out 
in the PK-3 Plus laboratory onboard the ISS under microgravity conditions. Details on 



the setup can be found in 18j. Dust particles injected into the main plasma formed a 
cloud around the center of the chamber. Several larger particles present in the chamber get 
sporadically accelerated and penetrate the cloud, thus generating projectiles [26|. Consider 



the experiment performed with argon gas (m n = 6.6 x 10~ 23 g). For the gas pressure 10 Pa 
and temperature T n = 300 K we have n n = 2.4 x 10 15 cm -3 and vt„ = 2.5 x 10 4 cm/s. 
Dust cloud was formed by melamine- formaldehyde particles with the diameter lad = 2.55 x 



10 4 cm and mass M d = 1.3 x 10 1 1 g 



19 



221 ]. therefore, v ~ 46 s 1 . The path of a 



subsonic projectile (u < c ~ 2.1 cm/s 19j, recorded with a high- resolution camera) revealed 
cavity deformation. The most obvious manifestation of the deformation is a visible shift of 
the projectile from cavity center to its front side. This shift can be accompanied by cavity 



10 

elongation (Fig. [2]). In the course of motion of the projectile, its velocity slowly increased 
from u ~ 0.7 cm/s to 1.4 cm/s. Cavity deformation becomes visible approximately in the 
middle of the path; the projectile velocity and acceleration at this point are 1.06 cm/s and 
2.6 cm/s 2 , respectively, the radius of a spherical cavity around the projectile at this point is 
R ~ 3.7 x 10~ 2 cm. 

One can crudely estimate the static pressure p st of dust cloud from the observed defor- 
mation threshold. According to the discussion at the end of Sec. HIl the projectile velocity at 
which a point of negative pressure emerges on the cavity surface, is lower than the velocity 
of deformation onset observed experimentally. Hence, in this case 2vR/9u CI ~ 1, and the 
system is at the boundary between the low- and high-velocity regimes. It is not surprising 
therefore that the corresponding estimates for the static pressure, K(= p s t/ n d) — M d uRu cr /2 
and K ~ 5Mdul T /8, both yield K ~ 3 eV, i.e., p st ~ 10 -7 Pa. Note that the correction for 
finite compressibility, u 2 R 2 /3c 2 , does not exceed 20%, and the correction for acceleration is 
negligible. 

Microgravity conditions can also be maintained on parabolic flights. Such experiments 
were performed using the IMPF-K2 chamber 
argon at 30 Pa and 2a d = 9.55 x 10" 4 cm (M d = 6.8 x 10~ 10 g and v ~ 33 s^ 1 ). Projectiles 
were accelerated by a special device up to supersonic or subsonic velocities. For a subsonic 
projectile, the cavity radius could be roughly estimated as R ~ 0.1 cm. The projectile 
velocity amounts to u ~ 1.8 cm/s, and cavity deformation is clearly seen in Fig. 3a of 



111 ] . The experiment was carried out with 



Ref. 



K on the size of a 



111 ] . Since we have no information concerning the dependence o: 
dust particle, we use the same value of K as for the experiment of Ref. 22| discussed above. 
Then the threshold velocity for this case does not exceed 1 cm/s. This corroborates the 
main notion of this study - the cavity around a moving projectile is normally deformed at 
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subsonic velocities. 

V. CONCLUSION 

We have investigated the possibility of deformation of a cavity around large projectiles 
slowly moving inside a cloud of small dust particles. This cavity is formed due to a strong 
Coulomb repulsion between the projectile and dust particles. Since the parameter (3d p char- 
acterizing this interaction is very large, cavity deformation occurs solely due to formation of 
a void adjacent to initially spherical cavity. This is similar to the process of void formation 
in a fluid under negative pressure. Since we approximate the dust cloud by a system of soft 
spheres, zero pressure is sufficient for void formation leading to cavity deformation. 

We model a collective particle subsonic motion in a complex plasma by nonviscous ir- 
rotational hydrodynamic motion of dust fluid and show that, in this model, the reason for 
cavity deformation is the friction between dust particles and neutral atoms of a carrier gas. 
It is this friction that stipulates emergence of a zone of negative pressure over the back of 
a cavity surface. If such zone covers substantial part of the surface, stall takes place, which 
is observed as cavity deformation. We have found that the deformation is more likely to be 
observed for larger projectiles. The deformation occurs when the projectile velocity exceeds 
some threshold subsonic velocity. 

Since typical projectile velocities are of the same order of magnitude as the sound ve- 
locity, we have studied the effect of finite compressibility of the dust fluid on velocity and 
pressure fields. With this purpose, we have obtained an exact solution of linearized gas dy- 
namics equations for compressible nonviscous fluid flowing about the cavity. In so doing, we 
included both the convective and friction terms in the Navier-Stokes equation. Calculations 
showed that corrections for finite compressibility are rather small. Generally, account of the 



12 

compressibility leads to a slight shift of the area covered by a negative pressure zone. 

Analysis of available experimental data has validated our approach. It was found that the 



cavity is deformed when either the velocity 22] or the size JjJ of a projectile are sufficiently 



large. It is worth mentioning that the phenomenon of cavity deformation is similar to the 
ormation of a void behind a cluster of smaller particles moving in the cloud of larger ones 



27| . For intermediate projectile sizes, a transition form spherical to deformed cavity was 
observed along the trajectory of an individual projectile. The transition threshold is in 
accordance with the developed theory. This made it possible to estimate the static pressure 
inside a dust cloud and, hence, opened up a possibility to use projectiles as a diagnostic tool 
for the dust equation of state. 

In this study, we restricted ourselves to investigation of the threshold of cavity deforma- 
tion and did not treat the shape of deformed cavity as well as displacement of the projectile 
away from the center of its cavity. This problem, besides all, would require an adequate 
account for local non-spherical charge distribution around the projectile. It is also possible 
that cavity deformation promotes the observed shift of the projectile from the center of a 



cavity, driven by an additional drag force exerted on the projectile 



23]. 



Appendix A: 

The solution of Eq. f[T2l) with boundary conditions ( 1131) is represented in the form 

Hx,Z) = i2fi(x)Pi(0, (Ai) 

1=0 

where fi{x) are the coefficients of expansion in the Legendre polynomials Pi{£,)- We substi- 
tute flSIl) in CL2D taking into account that V 2 P;(0 = ~K l + 1 ) x ~ 2p i(0 to reduce Eq. flTJJ) to 



13 



a set of ordinary differential equations in spherical coordinates for nonzero functions fi(x), 

H x2d £)- i{i+i)f,=q ' ix) ' (A2) 

where qo(x) = -kx~ 4 , q±(x) = ~^~ x ~ 5 ~~ ^ x ~ 8 ' 1^( x ) = ~ ra 1 + 2 Ka; ~ 4 ' ^ 3 ( x ) = ~^ x <2 + 



24 



-x 



-5 ° -8 



x . The boundary conditions for (1A2j) follow from (TT~3j) : 



5 2 

/l( °° ) = ' & 



0. (A3) 



The rhs of (IA2I) has the form J^bikX k , where are constants entering the functions qi(x). 

k 



Solution of (IA2I) is then given by 



fl{x) = ^ + ^ k(k-l)- 1(1 + 1) ^ (M) 



where C is defined by the second boundary condition (1A3|) . whence it follows that 

= E /(/ + 1) zib(5fc _ i) (tT!^ ~ i*) • (A5) 



Finally, we obtain 



6x 24x 4 
2 2 1 

A = + 



3x 2 ' 5x 5 12a; 8 ' 

f = k(- 1 I 1 
72 V6a; 6s 3 12x 4 

3 54 3 3 

J3 



(A6) 



■ 



5x 2 55x 4 5x 5 88a; 8 
Appendix B: 

The pressure correction due to finite compressibility can be found using the corresponding 
velocity correction, Eq. ( I14p . We rewrite Eq. (JSJ) in the form 

- c 2 V In p = i Vv 2 + vV<p, (Bl) 



14 

where <p = ip + ip\ and v = u + Vtp. By integrating Eq. (IBlj) and using Bernoulli's equation 
for compressible fluids we derive 

Po / vl \ po po fvl \ 2 

p ~ Po = — # ( y + V( pq j ~ ( v o ■ v i + v<pi) + ( y + V(po j ' ( B2 ) 

where po = p|r=oo and the terms up to (u/c) 4 are retained (for simplicity we assume 
c = const). We write the equation of state for the dust fluid in the form p — p st = c 2 (p — p ), 
which yields 

P° v o 1 ( Po^o . \ 2 mq x 

P = Pst ^Po<Po - Pov ■ vi - vp (pi + \ ~Y~ upoLp ° j ' ( J 

The first three terms of Eq. (1B3|) represent Eq. fl6]) with u = 0. Hence, by substituting in 
Eq. f TB3j) the expressions for <p and v [Eqs @ and (jSJ)] as well as for <pi and vi [Eq. ([Tj 
taken at the cavity surface, we get 

p u 4 



Pi(0 



.2 



-2 o-rf/rf 101 



c- L6 s V 440 

2025 / 4 5414 2 5237 \ 
1408 V ~ 337l> + 10125/ 



(B4) 
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FIGURE CAPTIONS 



Fig. [D Deformation of a spherical cavity (of the initial radius R) behind a projectile 
(large bullet in the center). The projectiles moves with the velocity u through a dust cloud 
(small bullets around) embedded in a carrier gas (gray background). 

Fig. [2J Almost spherical (a) and deformed (b) cavity around a projectile at different 
portions of its trajectory, observed with a high-resolution camera (inverted image). 
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FIG. 2. 



